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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  describes  the 
implementation  and  testing  of  waterborne  vessels  induced  pressure  and  stresses  into  the  three- 
dimensional  (3D)  shallow  water  module  (SW)  of  the  Adaptive  Hydraulics  (AdH)  numerical  code. 
A  description  of  the  inputs  required  by  the  user  is  also  provided  and  serves  as  the  user’s  guide. 

BACKGROUND:  The  movement  of  a  vessel  through  a  water  body  produces  currents  and  waves. 
Maynord  (1996)  provided  a  general  description  of  these  waves  and  currents  and  provided 
analytical  tools  to  describe  the  effects  of  these  currents  on  the  hydrodynamics  of  the  water  body. 
Stockstill  et  al.  (1999)  showed  that  by  applying  the  pressure  imposed  on  the  water  surface  by  the 
vessel  into  a  numerical  model,  these  currents  and  their  far-field  effects  can  be  accurately  simulated. 
In  this  technical  note  the  authors  expand  the  two-dimensional  (2D)  shallow  water  implementation 
presented  by  Stockstill  et  al.  (1999)  into  the  3D  shallow  water  realm. 

A  variety  of  3D  reservoir,  riverine,  coastal,  and  estuarine  simulation  models  exist  (such  as  ADH- 
SW3,  CH3D,  TABS-MDS,  etc.).  These  models  solve  some  form  of  the  Reynolds  averaged  Navier- 
Stokes  equations  using  the  hydrostatic  assumption.  The  governing  equations  for  velocity  (u,v,w) 
are  written  as 
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where,  Exx,  Exy,  Eyx,  Exz,  Eyy,  and  Eyz  are  kinematic  eddy  viscosities  (subscripts  indicates 
directionality),  po  is  the  reference  density,  and  P  is  the  hydrostatic  pressure.  The  horizontal 
( Exx ,  E  ,  E ,  E  )  and  vertical  ( Exz ,  Eyz )  eddy  viscosity  for  momentum  coefficients  are  utilized  to 
represent  the  Reynolds  stresses  and  are  obtained  through  turbulence  closure  schemes. 
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Under  the  hydrostatic  pressure  assumption,  vertical  accelerations  are  neglected,  and  the  equation 
for  w  yields  the  following  equation  for  the  pressure  term: 

v 

P(z)  =  Pa+f  9p{z)dz  (4) 

z 

where  Pa  is  the  atmospheric  pressure,  g  is  the  acceleration  due  to  gravity,  p{z)  is  the  density  at  a 
location  in  the  z  direction,  and  77  is  the  water  surface  elevation. 

INCLUSION  OF  VESSEL  PRESSURE  FIELD:  In  ADH-SW3,  the  pressure  applied  CP(z))  is  a 
function  of  the  atmospheric  pressure  (Pa),  usually  assumed  to  be  zero  in  the  absence  of  an 
atmospheric  pressure  field,  and  the  location  of  the  node  under  consideration  in  the  z  direction.  This 

hydrostatic  pressure  component  ( J  gp(zjdz  )  is  internally  normalized  by  the  product  of  gravity 

z 

and  density  so  that  the  resulting  pressure  is  in  terms  of  the  depth  of  the  node  under  consideration. 

Incorporation  of  boat  pressures  into  ADH-SW3  takes  advantage  of  this,  and  the  atmospheric 
pressure  ( Pa )  is  incremented  by  a  value  equal  to  the  draft  of  the  vessel  directly  over  the  node 
under  consideration.  This  results  in  Equation  4  being  modified  as  follows: 

P{z )  =  Pa+Pb+f  9P{z)dz  (5) 


Where  Pt  is  the  vessel  imposed  draft  or  pressure,  and  other  variables  are  the  same  as  defined 
previously. 

Figure  1  illustrates  this  vessel-imposed  pressure. 
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EXAMPLE:  The  incorporation  of  vessel  pressure  and  subsequent  stresses  was  tested  on  two 
problems:  (1)  a  rectangular  flume  oriented  at  a  45°  angle  to  the  x-axis  and  (2)  through  an 
application  to  Galveston  Bay. 

1.  Angle  flume.  The  first  test  case  is  a  rectangular  flume  with  dimensions  of  1000  m  in 
length,  50  m  in  width,  and  1  m  in  depth  (Figure  2).  The  vessel  has  a  draft  of  0.19  m  and 
other  properties  as  illustrated  in  Figure  3.  These  vessel  properties  are  defined  in  section 
“Boundary  File  Implementation.”  The  vessel  is  moving  through  the  flume  at  0.25  m/sec 
from  the  lower-left  to  the  upper-right  portion  of  the  flume. 


Figure  2.  Angle  flume  domain. 
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Figure  3.  Boat  definition  parameters  (parameters  are  defined  in  section  “Boundary  File 
Implementation”). 


Figure  4  illustrates  the  vessel-imposed  surface  pressure  field.  The  bottom  pressure  under  the 
vessel  footprint  should  be  1.19  m  because  of  the  additional  0.19  m  pressure  imposed  by  the 
vessel.  Figure  5  illustrates  the  bottom  pressure  under  the  vessel  footprint  and  shows  that  the 
model  correctly  computes  a  pressure  of  1.19  m. 
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Figure  4.  Vessel  pressure  imposed  on  the  surface. 


Figure  5.  Vessel  pressure  imposed  on  the  bottom. 


Figures  6  and  7  illustrate  the  effects  of  the  additional  pressure  imposed  by  the  vessel  on  the 
velocities  under  and  around  the  vessel.  It  is  seen  that  the  pressure  field  causes  a  velocity  response 
on  the  bed,  and  this  imparts  an  additional  stress  on  the  bed. 
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Figure  6.  Surface  velocity  response  to  vessel  pressure. 


Figure  7.  Bottom  velocity  response  to  vessel  pressure. 


2.  Galveston  Bay.  The  demonstration  of  boat  pressure  implementation  was  performed  on  a 
Galveston  Bay  example  (Figure  8).  This  demonstration  was  comprised  of  a  single  vessel 
drafting  1  m  traversing  the  Houston-Galveston  ship  channel  in  a  channel  of  approximately 
13  m  depth,  with  vessel  properties  shown  in  Figure  9.  These  vessel  properties  are  defined 
in  section  “Boundary  File  Implementation.”  This  example  was  run  with  no  tide,  wind, 
inflow,  or  other  forcings  other  than  the  moving  vessel.  Figures  10,  11,  and  12  illustrate  the 
vessel  pressure  field,  the  surface  velocity,  and  the  bottom  velocity,  respectively.  The  vessel 
is  moving  at  0.5  m/sec  into  Galveston  Bay,  from  lower  right  to  upper  left  of  the  figures. 
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Figure  8.  Galveston  Bay  ADH-SW3  domain. 


Figure  9.  Boat  definition  parameters  for  Galveston  Bay  vessel  (parameters  are  defined 
section  “Boundary  File  Implementation.”) 


Figure  10.  Vessel  pressure  footprint. 


6 


ERDC/CHL  CHETN-IV-106 
October  2015 


Figure  1 1 .  Surface  velocity  impact  of  vessel  pressure. 


Figure  12.  Bottom  velocity  impact  of  vessel  pressure. 


BOUNDARY  FILE  IMPLEMENTATION:  The  boat  pressure  field  is  activated  in  ADH-SW3  by 
the  inclusion  of  the  “OP  BT”  operation  parameter  card  in  the  ADH-SW3  boundary  conditions 
file  (*.bc).  In  addition,  a  boat  definition  file  (*.bt)  must  be  specified.  Presence  of  the  OP  BT  card 
in  the  boundary  conditions  file  and  the  absence  of  a  boat  definitions  file  will  cause  pre_adh  to 
exit  with  an  error. 

The  boat  definitions  files  have  a  series  of  cards  that  define  the  vessel  characteristics  such  as 
draft,  length,  etc.  The  first  line  in  this  file  indicates  how  many  vessels  are  being  simulated.  The 
card  descriptions  are  as  follows: 
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Boat  Definitions  Control  Cards 


BOAT 


Field 

1 

2 

FDEF 

Type 

char 

int 

Value 

BOAT 

>0 

BOAT  PARAMETERS 

Description 

Card  type 

Number  of  vessels  included 

SAILING  LINE  AND  INITIAL  VESSEL  POSITION 

Field 

Type 

Value 

Description 

1 

char 

FDEF 

Card  type 

2 

int 

>0 

Vessel  number 

3 

int 

>0 

Number  of  segments  in  vessel  path 

4 

real 

# 

x-coordinate  of  vessel’s  initial  position 

5 

real 

# 

y-coordinate  of  vessel’s  initial  position 

6 

real 

# 

Vessel  velocity  at  initial  position 

DRFT 

VESSEL  DRAFT 

Field 

Type 

Value 

Description 

1 

char 

DRFT 

Card  type 

2 

int 

>0 

Vessel  number 

3 

real 

# 

Vessel  draft,  units  of  length 

BLEN 

VESSEL  LENGTH 

Field 

Type 

Value 

Description 

1 

char 

BLEN 

Card  type 

2 

int 

>0 

Vessel  number 

3 

real 

# 

Vessel  length,  units  of  length 

BWID 

VESSEL  WIDTH 

Field 

Type 

Value 

Description 

1 

char 

BWID 

Card  type 

2 

int 

>0 

Vessel  number 

3 

real 

# 

Vessel  width,  units  of  length 

PBOW 

BOW  TO  LENGTH  RATIO 

Field 

Type 

Value 

Description 

1 

char 

PBOW 

Card  type 

2 

int 

>0 

Vessel  number 
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3 

real 

# 

Ratio  of  vessel  bow  length  to  vessel  length 

PSTR 

STREN  TO  LENGTH  RATIO 

Field 

Type 

Value 

Description 

1 

char 

PSTR 

Card  type 

2 

int 

>0 

Vessel  number 

3 

real 

# 

Ratio  of  vessel  stem  length  to  vessel  length 

CBOW 

FRACTION  OF  DRAFT  APPLIED  TO  PBOW 

Field 

Type 

Value 

Description 

1 

char 

CBOW 

Card  type 

2 

int 

>0 

Vessel  number 

3 

real 

# 

Ratio  of  vessel  bow  to  vessel  draft 

CSTR 

FRACTION  OF  DRAFT  APPLIED  TO  PSTR 

Field 

Type 

Value 

Description 

1 

char 

CSTR 

Card  type 

2 

int 

>0 

Vessel  number 

3 

real 

# 

Ratio  of  vessel  stem  to  vessel  draft 

PROP 

PROPELLER  PARAMETERS 

Field 

Type 

Value 

Description 

1 

char 

PROP 

Card  type 

2 

int 

>0 

Vessel  number 

3 

int 

1  or  2 

Propeller  type 

4 

real 

>0 

Propeller  diameter 

5 

real 

>0 

Distance  between  centers  of  propellers 

6 

real 

>0 

Tow  vessel  length 

7 

real 

>0 

Distance  from  props  to  stem  of  tow  vessel 

SDEF 

SAILING  LINE  SEGMENT  DEFINITION 

Field 

Type 

Value 

Description 

1 

char 

SDEF 

Card  type 

2 

int 

>0 

Vessel  number 

3 

int 

>0 

Sailing  line  segment  number 

4 

int 

0  or  1 

Type  of  segment  (0=line;  l=arc) 

5 

real 

# 

x-coordinate  of  segment  end 

6 

real 

# 

y-coordinate  of  segment  end 

7 

real 

# 

Vessel  velocity  at  segment  end 
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8 

9 

10 

ENDD 

real 

real 

real 

# 

# 

+1/-1 

x-coordinate  of  arc  center 
y-coordinate  of  arc  center 

Direction  of  arc  turn  (+1  for  left,  -1  for  right) 

BOAT  PARAMETERS 

Field 

Type 

Value 

Description 

1 

char 

ENDD  0  0 

Card  type 

SUMMARY  AND  CONCLUSIONS:  This  CHETN  describes  the  development  and  testing  of  the 
capability  to  simulate  vessel-generated  currents  and  stresses  into  ADH-SW3.  Two  examples  are 
provided  demonstrating  this  capability,  and  these  examples  demonstrate  that  the  ADH-SW3  code 
is  accurately  simulating  the  resultant  pressures  and  currents.  This  implementation  of  vessel¬ 
generated  currents  and  stresses  ignores  the  vertical  accelerations;  therefore,  velocities  near  the 
vessel  are  not  accurate  and  should  not  be  used  for  the  design  of  engineered  channels  without 
careful  engineering  judgment.  The  “Boundary  Condition  File  Implementation”  section  also 
provides  all  the  information  necessary  for  a  user  to  include  and  activate  this  capability  in  ADH- 
SW3. 

POINTS  OF  CONTACT:  This  CHETN  was  prepared  as  part  of  the  Flood  and  Coastal  (F&C) 
Damage  Reduction  program  of  the  U.S.  Army  Corps  of  Engineers.  The  POC  for  technical  inquiries 
is  Jennifer  N.  Tate  (Jennifer. N.  Tatefa), USA CE.Armv.mil).  For  information  about  the  F&C  program, 
please  contact  the  F&C  Program  Manager,  Cary  A.  Talbot  (CaryA.  Talbotfw, usace.armv.mil).  This 
technical  note  should  be  referenced  as  follows: 

Savant,  G.,  and  C.  J.  Trahan.  2015.  Shallow  Water  Adaptive  Hydraulics  (ADH- 
SW3):  Waterborne  vessels.  ERDC/CHL  CHETN-IV-106.  Vicksburg,  MS:  U.S. 

Army  Engineer  Research  and  Development  Center. 

An  electronic  copy  of  this  CHETN  is  available  from  http: //chi,  erdc.  usace. annv.mil/chetn . 
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